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ABSTRACT 

Context. Molecular outflows from very low-mass stars (VLMSs) and brown dwarfs have been studied very little. So far, only a few 

CO outflows have been observed, allowing us to map the immediate circumstellar environment. 

Aims. We present the first spatially resolved H 2 emission around IRS54 (YLW52), a ~0. 1-0.2 M Q Class I source. 

Methods. By means of VLT SINFONI K-band observations, we probed the IF emission down to the first ~50 AU from the source. 

Results. The molecular emission shows a complex structure delineating a large outflow cavity and an asymmetric molecular jet. 

Thanks to the detection of several IF transitions, we are able to estimate average values along the jet-like structure (from source 

position to knot D) of A v ~28mag, T ~2000-3000K, and IF column density N(H 2 )~1.7xl0 17 cirr 2 . This allows us to estimate a 

mass loss rate of ~2x 10~ 10 M G yr~' for the warm IF component . In addition, from the total flux of the Br y line, we infer an accretion 

luminosity and mass accretion rate of 0.64 F Q and ~3xl0~ 7 M Q yr -1 , respectively. The outflow structure is similar to those found 

in low-mass Class I and CTTS. However, the F act ./Ff, D ; ratio is very high (~80%), and the mass accretion rate is about one order of 

magnitude higher when compared to objects of roughly the same mass, pointing to the young nature of the investigated source. 



Key words, stars: formation - stars: circumstellar matter - 
IRS54- Infrared: ISM 



1. Introduction 



ISM: jets and outflows - ISM: individual objects: YLW52, ISO-Oph 182, 



> 

o 



x 



Protostellar jets and outflows are associated with the first stages 
of stellar evolution. They are usually found in active young stel- 
lar objects (YSOs), in which a significant fraction of the cir- 
cumstellar material is still accreting onto the protostar. Jets and 
outflows can therefore be considered an outcome of the accre- 
tion activity. During the early evolutionary phase, outflows are 
usually traced through molecular emission lines. Molecular out- 
flows from low- and high-mass protostars ha ve been extensively 
studie d through observa tions of CO lines (Bachille r & Tafallal 
1999; lArce et al.1 12007b . These observations have revealed 
many common observational characteristics between both mass 
regimes, which points to a commo n outflow mechan ism working 
from low- to high-mass protostars (lArce et al . 20071). 

Molecular outflows from VLMSs and brown-dwarfs (BDs) 
have been studied very little, especially during their early evo- 
lutionary phase. Only recently have CO submillimetre obser- 
vations been able to provide direct imag ing of a small num- 
ber o f molecular outflows from VLMSs (iPhan-Bao et al.ll2008l 
1201 lb . Most of the few observations of jets from BDs and 
VLMSs, however, involve relatively evolved, classical TTauri- 
like YSOs, mainly studied through detecting forbidden emis- 
sion lines (FELs) dFernandez & Comeronl l2Q0lt IWhelan et al.l 
1201 2|) in their spe ctra, and/or spectro-astrometry of FELs 
(IWhelan et all [2005. and ref. therein). So far, no observation of 
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resolved molecular hydrogen emission line (MHEL) regions has 
existed for VLMSs or BDs. The presence of MHEL regions are 
typical of the spectrum of low-mass Class I source s dDavis et al.l 
2001: iGarcia Lopez et alJl2008t iDavis et alJl201 lb . and they are 
usually associated with an FEL region. Thus, one should also 
expect MHEL regions around young BDs and VLMSs. 

In this context, we present here the first H2 1-0S(1) spectro- 
imaging of an outflow from a Class I VLMS, IRS54 (YLW 52). 
This source (a=16:27:51.7, £=-24:3 1:46.0) is located outside of 
the main clouds in the Ophiuchus star-forming region, and it has 
been classified as a late-sta ge Class I source with a bolomet- 
ric luminosity of -0.78 L G (Ivan Kempen et alJl2009l) . IRS54 is 
thus in its main accretion phase, representing one of the lowest 
luminosity sources for which an H2 outflow has been spatially 
resolved. 



2. Observations and data analysis 

The data were acquired on 14 June 2010 at the Very Large 
Telescope at Paranal Observatory, Chile, using the integral field 
spectrograph SINFONI at medium resolution in the K-band 
(R~4000). The chosen pixel scale was 100 mas, corresponding 
to a 3"x3" field of view. The observations were acquired under 
0'.'4 seeing (DIMM FWHM), leading to a spatial resolution of 
48 AU at the location of IRS54 (d~120pc). The total integration 
time is 1200 s. To correct for atmospheric response, observations 
of a telluric standard star of spectral type B were performed. 
The main data reduction process was done using the SINFONI 
data-reduction pipeline, i.e., dark and bad pixel masks, flat-field 
corrections, optical distortion correction, and wavelength cali- 
bration using arc lamps. To test the goodness of the wavelength 
calibration, we applied the final wavelength transformation ma- 
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Fig. 1. Continuum-subtracted H2 1-0S(1) images of IRS54. Panels a 
and c: average over two spectral channels at -64kms~' and -30kms~' 
(panel a), and +40kms~' and +75kms~' (panel c). Panel b: Single 
spectral channel at +6kms~'. The velocities are corrected from an av- 
erage cloud velocity of 3.5 kms~' dWouterloot e t al. 2005; Andr e et all 
2007). Panel d: average H 2 1-0S(1) emission over the previous veloc- 
ity channels. Overplotted are the locations of six spatial regions where 
spectra were extracted. For reference, contours of the continuum (near 
the H2 2.122 fim line) down to the FWHM size have been overplotted at 
the centre of every image (dashed-blue contours). Contours in all panels 
show values of 1.3, 2.9, 4.5, 5.3, 6.1, 7.7, 9.3xl0~ 14 erg s _1 cnT 2 /rni~'. 



trix to a sky cube. By measuring the wavelength of the OH lines 
present in the cube, we found a systematic wavelength shift of 
~2.2A with respect to their theoretical value. The final data cube 
was thus shifted in wavelength to take this error into account. 
In addition to the SINFONI pipeline, the STARLINK software 
was used to correct the spectrum from atmospheric absorption 
and to flux-calibrate the data. With this aim, the standard star 
spectrum was extracted by collapsing the central region of the 
data cube. The Bry line was then removed from the standard 
spectrum before dividing it by a normalised blackbody at the ap- 
propriate temperature. Then, the spectrum was grown up into a 
60x70 pixel cube, as was that of our SINFONI data, to correct 
for the telluric features and to flux-calibrate the data. 

The brightest lines present in the K-band SINFONI spectrum 
are H2 lines from different rovibrational levels and the Br y line. 
Because both continuum and line emissions are observed simul- 
taneously, an accurate subtraction of the continuum can be per- 
formed, allowing us to obtain very clean images of the different 
line-emitting regions. To better probe the jet morphology close 
to the central source, the continuum emission was removed from 
the IFS cube using the IRAF subroutine "CONTINUUM" iter- 
atively along the dispersion axis at each spatial position of the 
datacube. All images shown in this letter have been continuum- 
subtracted using this method. In addition, integral field spec- 
troscopy enables us to obtain precise positional measurements 
of emission knots relative to each other and to the source con- 
tinuum. For instance, a ~ 16 mas positional accuracy can be ob- 
tained by taking our seeing conditions and a signal-to-noise ratio 
of only ~ 10 into a ccount (position accuracy ~Seeing/(2.35xSN); 
iDaviset alJl20TTh . 

3. Outflow physics and morphology 

Integral field spectroscopic observations allow us to retrieve di- 
rect information about the morphology and kinematics of the 
H2 emission close to the Class I protostar IRS54. Indeed, the 
H2 2. 122 fim line is one of the brightest lines tracing outflow ac- 
tivity in the K-band. 

FigureQ] shows the averaged continuum-subtracted H2 1- 
0S(1) emission across the datacube (panel d), along with the 
blue- and red-shifted H2 outflow components (panels a and c). 
In addition, the H2 emission at rest velocity (corresponding to 
one pixel in the datacube dispersion direction) is shown in panel 
b. The H2 1-0S(1) spectral image in panel d was contructed by 
averaging five different velocity channels at the local standard of 
rest (from — 64 km s -1 to +75kms _I ), while the blue- and red- 
shifted images correspond to the averaged emission over two 
pixels in the dispersion direction (-64km s _1 and -30km s _1 , and 
at +40kms _I and +75kms~ I ). The figure shows a very com- 
plex H2 morphology with gas displaying an X-shaped spatial 
distribution superimposed on a more collimated structure, possi- 
bly jet-like (see discussion bellow), located westward of IRS54 
and extending over ~1" (~120 AU). All regions show a knotty 
structure, with condensations A, B, C, E, and D displaced from 
the source (Aa,AS) at about (l'.'3,-0'.'4), (0'.'4, 075), (073, 074), 
(-071,0), and (-073,0), respectively (see panel d). The overall 
gas distribution is very similar to the one observed in molecu- 
lar jets from embedded low-mass p rotostars, as usuall y traced 
by CO/molecular observations (e.g.. iTafalla et al.ll2004l) . In this 
case, however, two well-defined red- and blue-shifted lobes can- 
not be identified. As shown from the velocity channels repre- 
sented from panels a to d, red- and blue-shifted emission can be 
roughly associated with the same spatial regions. It should be 
noted that the velocities of the channel maps do not correspond 
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to the H2 2.122yum line peak velocity at any spatial position. In- 
deed, this line shows a very broad profile with a full width zero 
intensity of ~200kms~ 1 , but with a peak velocity of ~0kms _1 
(see Fig. 0. This may suggest that the outflow lies very close to 
the plane of the sky. 
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Fig. 2. Line profile of the H2 2.122yum line extracted for six diffe- 
rent spatial positions (see Fig.[T]and Table[T]for information about the 
spatial distribution of the different extracted areas). All line profiles are 
normalised to the flux of knot D. 

In addition, thanks to detecting several H2 lines, excita- 
tion diagrams could be constructed to derive the temperature, 
H? column density, and ex t inction of the molecular gas (e.g., 
Cara tti o Garatti et alj |2006; Garc ia Lopez et al1l2010l) . The re- 
sults are shown in TableQ] The H2 line fluxes (see Table lA. It 
were measured at each of the positions indicated in Fig. [T] by 
simulating "long-slit" spectra with a fixed slit width correspond- 
ing to an extraction aperture of 3x3 pixels. Temperature and 
N(H 2 ) values range from -1840K to 3340 K and from 2.4 to 
20xlO I6 cirT 2 . The visual extinction, temperature, and NCH2) 
are higher westward of IRS54, increasing along the jet-like 
structure towards the source position (knot D and on source), 
which suggests that the western lobe is the red-shifted lobe of 
the outflow and that the X-shaped and jet-like structures have di- 
fferent physical conditions. The extinction values found towards 
the source are very similar to those rep orted in previous w orks, 
where Ay values >25 mag are reported (Hai sch et al.l l2004). 

Finally, we note that the Br 7 emission is not spatially re- 
solved, and it matches the stellar continuum. 

4. Accretion and ejection properties 

From our observations, it is also possible to derive the ac- 
cretion luminosity and the mass accretion and ejection rates 
of IRS54. Using the relation between the luminosity of the 
Bry lin e (L(Bry)) and the accretion luminosity (L acc ) de- 
rived by ICalvet et aD (|2004), we found an accretion luminos- 
ity of L flcc ~O.64L . The luminosity of the Bry line was mea- 
sured from the integrated flux across the Br y line in our cube 
(F=4.08xl0~ 14 ergs~' crrr 2 ) and corrected by a visual extinc- 
tion of 30 mag. Consideri ng the modest bolometri c luminosity of 
this source (Li o ;~0.78Lo;|van Kempen et ail2 009). the accretion 
luminosity is very high, with an L flcc /L/, /value of -80%, which 
is consistent with a very young and a ctive YSO. This supports 
the idea that IRS54 is a young YSO (Ivan Kempen et alJ l2009) 



Table 1. H 2 diagnostics. 
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still accreting large amounts of circumstellar matter. Assum- 
ing L bo i 

— L fl cc+L*, we find a stellar luminosity of L* — 0.14Lq. 
We derive a stellar mass and radius of 0. 1-0.2 M and ~2R by 
placing the stellar luminosity into an HR diagram an d consider- 
ing the bit hline locus usin g the evolutionary tracks of Siess et all 
(120001) and lStahlerl (119881) . This indicates that IRS54 has spectral 
type M and is a very low-mass star. 

A mass accretion rate of ~3.OxlO _7 M yr _1 is in- 
ferred from the observed L acc , following the expression 
Mfl CC =(L acc R./GM H ,)x(l-R./R i -), R , being the inner radi us of 
the accretion disk (with R,=5 R«, see lGullbring et al.lll998h . The 
derived value is higher than those found in more evolved object s 
of roughly the same mass dNatta et al.ll2006t iGatti et al] |2006). 
again pointing to the young nature of this source. 

The mass-loss rate carried by the warm molecular com- 
ponent (Mh 2 ) can t> e computed from the derived NTHz) val- 
ues using the exp r ession Mh 2 -2 H m#N(H2) A dv,/dl r (see, e.g., 
iDavis et al]l200ll 120111) . Here, A is the area of the emitting 
region, fi is the mean atomic weight, and dl, and dv, are the 
projected length and the tangential velocity. We have assumed 
that the area where the molecular hydrogen is emitted is equal 
to the extent of the flow along the jet axis from position S to 
D (~0'.'5) multiplied by the width of the flow (i.e. the seeing). 
We have considered an average N(H2) value for this region of 
~1.7xl0 17 cm~ 2 . To retrieve y t , the inclination angle of the out- 
flow with respect to the line of sight should be known. As out- 
lined before, the radial velocities suggest that the outflow is al- 
most in the plane of the sky, and thus, an estimate of the jet ve- 
locity from an unknown flow angle is too uncertain. Therefore, 
we assume the width of the H2 lines presented in Fig.|2]as a lower 
limit to the jet velocity, i.e. ~200kms~'. This gives a value of 
> 1.6x1 0~ 10 M yr~' which is about two orders of magnitude 
lower than expected if compared with the derived M acc value 
(M ol i,/M acc ~Q.l). The computed value is also around two to 
three orders of magnitude lower than the M# 2 value derived for 
low-m ass Class I sources using the same technique (IDavis et al.l 
|20TT1) . However, it is worth noting that these sources are likely 
to be more massive than ours. This suggests that (1) most of 
the outflow material is transported by a cooler and denser com- 
ponent than trac ed by the near-IR H2 lines, such as the H2 pure 
rotational lines ( Caratti o Garatti et al . 2008); or (2) most of the 
H2 has been dissociated, and the jet is mainly atomic. The latter 
has been observed in other Class I jets, where the measured Mh 2 
is around one or der of magnitude l ower than measur ed from the 
[Fen] emission jDavis et alj|201 lhlNisini et alj|2005h . 



5. Discussion 

Molecular outflows, H2 jets in particular, have been observed 
from Class to Class II sources through a wide range of masses. 
Here, we show that MHEL regions are also common in the first 
evolutionary stages of VLMSs, suggesting that the same launch- 
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ing mechanism is acting from very low- to high-mass stars and 
across a wide range of evolutionary stages. 

In general, MHEL properties are consistent with shock- 
heated gas from the inner regions of Herbi g-Haro objects or 



from spatially extended wide-angle winds (Davis et al 



12001b 

ICaratti o Garatti et al.l2006tlBeck et al.l2008l;lDavis et al.l201 ll) 
The reported gas temperatures and excitation properties of 
IRS54 are consistent with shock-heated material, as well. Simi- 
lar valu es have also been found in low-mas s Class I sources and 
CTTSs dBeck et al]2008l:lDavis et al.l201 lh . The overall outflow 
structure can be interpreted as the sum of H2 emission excited 
along a wide-angle cavity, plus the contribution from a jet-like 
structure. This latter structure is, however, only observed in one 
of the lobes. This probably indicates different excitation condi- 
tions for the two lobes of the jet, i.e. an asymmetric jet in which 
the blue-shifted lobe might have a higher velocity, dissociating 
the molecular H2 component. If this is the case, ionic/atomic 
emission should be present in the eastern jet. Velocity asymme- 
tries have bee n observed before in o ther protostellar jets, such as 
RW Aur (e.g JMelnikov et ai1l2009l) . 

IRS 54 is known to possess a large-scale S-shaped bipolar 
jet as traced by a chain of H2 knotfl extending over several 
arcminutes (feature f08-01 in Khanzadvan et al.l 120041) . From 
wide-field H2 images, the jet appears to precess with a position 
angle (PA) changing from ~72° in the outermost knots passing 
through ~78° in the middle part of the jet dow n to 90° closest 
to the source (MHO 2 129; see also Fig. B.l in J0rg ensen et al.l 
2009). The strong jet precession may then be at the origin of the 
outflow cavity, suggesting that the low-velocity H2 emission is 
gas excited along the cavity walls by the interaction of a wide- 
angle wind with the ambient material. On the other hand, the 
PA of the internal jet-like structure detected in our H2 spectral 
images is consistent with the PA of the innermost knots of the 
large-scale H2 jet (i.e., MH02129). This may support the idea 
that the denser and more reddened jet-like structure in our im- 
ages is indeed the upstream region of the H2 jet observed in the 
wide-field images. Jet precession alone can, however, hardly ac- 
count for all the observed features shown in Fig.Q] By assuming 
a lower limit to the jet velocity of ~200kms _1 , the dynamical 
age of condensations A, B, C, D, and E are around 3.9, 1.9, 1.4, 
0.9, and 0.3 yr. This would lead to a precession period of about 
four years in order to generate the structures (A, B, and C) at 
each side of the so-called cavity (opening angle ~65°). This pe- 
riod is t oo short when compared with typical jet fast prece ssion 
periods dRosen & Smithll2004t ICaratti o Garatti et al.ll2 008). and 
not consistent with the longer period suggested by the l arge-scale 
jet images (see f08-01 field in lKhanzadvan et alJl2004l) . 

On the other hand, a wide-angled wind might be able to ex- 
plain both the broad line profiles presented in Fig.|2]and the ob- 
served features in Fig.Q] Considering the jet-like structure alone 
(condensations E and D), we have measured the H2 jet width as 
a function of the distance from the source (see, Fig. IB. 11 . We es- 
timate a full opening angle of the flow of ~23 °, fittin g a straight 
line to the points from 40 to 140 AU in Fig. lB.lK see lDavis et al.l 
l20Tll for more details). Similar values are also found in jets 
from low-mass Class I sources and CTTSs, which show opening 
angles between 20° and 42° ( Hartiganet al.1 l2004t [Davis et al. 



1201 li) . Beyond >50 A U, the jet width slow l y increases with d is- 
tance (see Fig. lB.ll and lHartigan et alj2004tlDavis et alj201 lh i 



expec ted for a free lateral expansion of a supersonic jet (jCabrit 
120071) . Therefore, the jet collimation must take place within the 
first <50AU from the source, in agreement with MHD wind 



named MHO 2128-2131, see http://www.astro.ljmu.ac.uk/MHCat/ 



models dDougados et alJ2004h . iPanoglou etall d2012h show that 
molecular hydrogen can survive along MHD disk-wind stream- 
lines. The presence of H2 emission detected at only ~48 AU 
from the source, together with the broad line profiles shown in 
Fig.|2j tend to favour the presence of an MHD disk-wind model. 
Nevertheless, a wide-angled wind alone cannot account for the 
precessing jet observed in the large-scale images. Therefore, a 
combination of both a wide-angled wind and a precessing jet is 
possibly the main mechanism behind the complex H2 structure 
in IRS54 (Fig.EIl). 

6. Conclusions 

In this letter, we present the first spatially resolved H2 emission 
(MHEL) region around IRS54, a Class I VLMS. The H 2 emis- 
sion was detected down to the first ~50 AU from the source, and 
it shows a very complex morphology. The emission might be 
interpreted as coming from the interaction of a wide-angle wind 
with an outflow cavity and a molecular jet. In addition, the de- 
tection of several H2 line transitions and the Br y line allows us 
to derive various accretion/ejection properties: 

- We computed the extinction, H2 column density, and tem- 
perature values at different outflow spatial positions from the 
analysis of excitation diagrams. The highest values are found 
westward of IRS54, increasing along the jet-like structure to- 
wards the source position. Average values of A„ ~28 mag, 
T=2000-3000K, and N(H 2 )~1.7xl0 17 cirT 2 are found along 
the jet-like structure (knots D and S). 

- We inferred an accretion luminosity and mass accretion rate 
of 0.64 L and 3xl0" 7 M Q yr" 1 from the total flux of the Bry 
emission. From the computed L acc and the L/, c / value found 
in literature, we derived L» ~O.14L , which yields a stellar 
mass of ~O.l-O.2M .The accretion luminosity accounts for 
~80% of the total luminosity. This, together with the high 
M acc value, points to the young nature of IRS54. 

Acknowledgements. We thank the anonymous referee for the comments that 
helped improve the paper. 
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Appendix A: H 2 line fluxes and spectra 
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Fig. A.l. Continuum-subtracted spectra at five different positions along the outflow (see Fig.[T]and Tab lA. 1 \ . 



Table A.l. H 2 line fluxes. 
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Appendix B: Outflow width and morphology 
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Fig. B.l. Width of the H 2 1-0S(1) jet-like structure as a function 
of the distance from the source. The jet width has been estimated 
by extracting 2-pixel-wide vertical slices to the average H 2 spectral 
image presented in Fig.[T]and fitting a single Gaussian function. Zero 
offsets corresponds to the source positions. 




Fig. B.2. Sketch showing IRS54 outflow/jet morphology. 



